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PHOTOPRODUCTION OF HEAVY QUARKS AT HERA 

Talk given at the 3 rd UK Phenomenology Workshop on HERA Physics, Durham (UK), 
September 2(f h -2$ h , 1998. 

Heavy-flavoured hadrons are mainly produced at HERA in those events in which 
the virtuality of the photon exchanged between the electron and the proton is very 
small (see R. Thome's talk [0 for a review on heavy quark production in deep-inelastic 
scattering.) The electron, in this case, is equivalent to a beam of on-shell photons 
carrying a fraction of the electron energy distributed according to the Weizsacker- 
Williams function p|. The underlying production mechanism is therefore a photo- 
production one, which has been studied extensively in fixed-target experiments (see 
ref. || for a complete bibliography); HERA offers the opportunity of testing the pre- 
■^j- . dictions of perturbative QCD for this process in an energy regime which is totally 

unexplored, the available centre-of-mass energy being about one order of magnitude 
larger than in fixed-target experiments. This is an interesting opportunity, since pho- 
toproduction is much less affected by nonperturbative and higher order effects than 
hadroproduction. 

A well-known feature of photoproduction processes is that the incoming photon 
may fluctuate into a hadronic state before undergoing a hard collision. The corre- 
sponding contribution to the cross section is referred to as hadronic (or resolved), to 
contrast with the case when the photon directly interacts with the hadron [point-like 
or direct.) Therefore, a differential photon-hadron cross section can be written as the 
sum of a point-like and a hadronic photon contribution: 

b : 

da^\P„ P H ) = da^ t (P 7 , P H ) + da^(P„ P H ) . (1) 
According to factorization theorems, the two contributions can be expressed as 

d Vpofl( P i, p b) = J dxfj"' \x, n F )da lj {P 1 ,xP H ,a s {n 

), (2) 

j 

x da ij {xP 1 , yP H , a s (/4)> /4> a4> /%) > ( 3 ) 



Oh 



where f\ and /j are the parton densities in the photon and in the hadron re- 
spectively. The point-like and the hadronic components of the photoproduction cross 
sections are very closely related, and only their sum is physically meaningful. In fact, 
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Figure 1: Total cross section for the photoproduction of cc pairs, as a func- 
tion of the "fp centre- of-mass energy: next-to-leading order QCD predictions 
versus experimental results. 



the separation of a cross section into a point-like and a hadronic component is am- 
biguous beyond leading order: different choices of the factorization schemes lead to 
different definitions of the two components. This is explicitly shown in eqs. (|2|) and 
(|3|): the same arbitrary scale /x 7 , related to the subtraction of collinear singularities 
of quarks emitted by the incoming photon, appears in both equations. 

The parton densities in the photon are quite soft. Therefore, the hadronic com- 
ponent is only important for large CM energies and small masses for the produced 
system, which is precisely the case for charm production at HERA. This process can 
therefore be used to constrain the densities in the photon, which are experimentally 
very poorly known at present. 

In fig. |I] I present measured values of the total charm cross section in fixed tar- 
get and HERA experiments, superimposed on the theoretical predictions, for three 
different choices of the parametrization of the photon parton densities. Each pair of 
curves represents the change in the theoretical predictions due to variations of the 
renormalization and factorization scales; the corresponding uncertainty is quite large, 
due to the small value of the charm quark mass, which sets the energy scale of the pro- 
cess. It is also clear that the impact of the hadronic component is negligible at fixed 
target energies, while it becomes important in the HERA energy range. The data 
show a satisfactory agreement with the curves obtained by choosing m c = 1.5 GeV 
and MRSG Q and GRV-HO || for the partonic densities in the proton and the pho- 
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ton respectively. Remarkably, a single choice of the input parameters allows a good 
description of the data in the whole energy range considered. 

It is apparent from fig. [I] that no definite conclusion can be drawn on the photon 
densities. This is partly because of the presence of inconsistencies in the low-energy 
data, and partly because of the intrinsic uncertainties of the HERA measurements. In 
fact, the experiments are only sensitive to production in the central rapidity region; 
furthermore, a small-p T cut is applied to the data, in order to clearly separate the 
signal from the background. This kinematical region typically involves small- x values, 
and the extrapolation is therefore subject to uncertainties due to our ignorance of the 
small- a; behaviour of the parton densities. 

Single- inclusive distributions have also been measured by ZEUS and HI; a detailed 
description of the data compared to theoretical predictions is given in the contribution 
of J. Cole to these Proceedings ||. A general comment to these studies is that 
the agreement appears to be satisfactory; however, it is necessary to wait for more 
statistics to draw significant conclusions. In fact, if data with larger statistics were 
available, it would be possible to consider the production processes initiated by very 
energetic photons only, instead of integrating over the whole photon spectrum using 
the Weizsacker- Williams distribution. 

At HERA energies, large terms proportional to powers of log(S/m^) appear in 
the charm cross section, that may spoil the convergence of the perturbative series. 
The problem of resumming these terms (small-x effects) has been studied in ref. [0 
in different contexts. In ref. || it has been estimated that these effects are negligible 
with respect to the other sources of uncertainty on the theoretical predictions. 

A similar problem arises when the transverse-momentum distribution is consid- 
ered; in this case, terms proportional to powers of \og(p T /m c ) are present in the 
coefficients of the perturbative expansion, which are potentially dangerous at large 
p T . These logarithms can be resummed by observing that, at high p T , the heavy-quark 
mass is negligible, and by using perturbative fragmentation functions H [10], 11]. It 



has been checked that the fixed-order and the resummed results of |J agree in a very 
wide range in p T . Recently, a formalism has been developed [[l!| that allows matching 
the two approaches (the fixed-order calculation with m/0 and the resummed result 
for massless quarks). It has however been applied only to hadroproduction processes. 

Let us now briefly consider b production. This case is more reliably described in 
perturbative QCD, because of the larger value of the 6-quark mass. All the theoretical 
uncertainties are strongly reduced with respect to the charm case. Also, the hadronic 
component of the photon is less important in b production, which involves larger values 
of x. It should be remarked that a statistically significant study of bottom production 
will be possible at HERA only with a luminosity upgrade (the total rates are about 
a factor 200 smaller than those of charm production), and in any case a comparison 
with theoretical predictions will only be possible considering electroproduction in the 
Weizsacker-Williams approximation. 

Figure shows that the transverse momentum distribution of b quarks at HERA 
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Figure 2: Full uncertainty on the transverse-momentum distribution for 
bottom electroproduction (Weizsacker- Williams approximation) with Peterson 
fragmentation and a pseudorapidity cut. 



can be predicted by perturbative QCD quite accurately. Even with the LAC1 13 
set the hadronic component affects the prediction only marginally; this is basically a 
consequence of the applied pseudorapidity cut. 

The first measurements of the b production cross section at HERA have recently 
been presented by the HI Collaboration Hl4|. 

If the integrated luminosity of the HERA collider will increase to 100 pb _1 or more, 
it is conceivable to perform analyses of double-tagged charm events |15|], and therefore 



to study double-differential cross sections. This would be a precious opportunity in 
many respects; for example, a direct determination of the proton gluon density f^ 
16| , [17| would become possible. 



In photon-hadron or electron-hadron collisions, the gluon density enters the cross 
section in a simpler way than in the case of hadronic collisions. If one could determine 
the invariant mass and rapidity of the produced system, the leading-order cross section 
would be directly proportional to f^ and to other calculable factors. For example, the 
LO heavy-flavour electroproduction cross section can be written, in the Weizsacker- 
Williams approximation, and neglecting the hadronic component, in the following 
way: 



x 



i- &\x y ,4)&\x g ,MKM^), (4) 



dy^dM^ 9 dx 9 dM^ E* J i WW wwjg^J, 

where M Q q is the invariant mass of the heavy-quark pair, and y Q Q is the rapidity of 
the pair in the electron-proton centre-of-mass frame (positive rapidities are chosen in 
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the direction of the incoming proton); E = \/S is the electron-proton centre-of-mass 
energy, and 



exp(-?/QQ), 



E 



M. 



E 



exp(y Q g). 



(5) 
(6) 



The function is the Weizsacker- Williams distribution. All quantities on the right- 
hand side of eq. (||) are calculable, except for fjf \ which can therefore be measured. 
The inclusion of radiative corrections up to NLO is straightforward, as explained in 



M c - 5 >15 GeV 



Solid: pointlike 
Dashed: hadr. (GRV)j 
Dotted: hadr. (LAC1) - 





Figure 3: x g distribution in ep collisions (Weizsacker-Williams approxima- 
tion) at HERA, for m c = 1.5 GeV and the MRSG parton densities in the 
proton. 



ref. [0. 

The method outlined above for the extraction of the gluon density in the proton 
will only work in those kinematical regions where the hadronic component is sup- 
pressed. This possibility is illustrated in fig. |3|, where the NLO QCD distribution 
in the variable x g in electron-proton collisions at = 314 GeV is presented. The 
partonic densities in the proton are given by the MRSG set, while both the LAC1 
and GRV-HO sets for the photon are considered. Figures ^a)-0c) show the effect of 
applying a cut on the invariant mass of the pair. Even in the case of the smallest 
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invariant-mass cut, there is a region of small x g where the hadronic component is 
actually negligible. Increasing the invariant-mass cut, the hadronic component de- 
creases faster than the point-like one. This is because for large invariant masses the 
production process of the hadronic component is suppressed by the small value of 
the gluon density of the photon at large x. By pushing the invariant-mass cut up to 
20 GeV, it turns out that the point-like component dominates over the hadronic one 
for x g values as large as 10 _1 . Therefore, it can be concluded that the theoretical 
uncertainties affecting the charm cross section, in the range of 10 -3 < x g < 10 _1 , are 
small enough to allow for a determination of the gluon density in the proton by using 
invariant-mass cuts to suppress the hadronic component. Figure |3]d) shows the effect, 
on the x g distribution, of a small-p T and a pseudorapidity cut applied to both the 
charm and the anticharm, as usually done in HERA analyses. It is interesting to note 
that in this case the point-like component is dominant in the whole kinematically 
accessible range, even without applying an invariant-mass cut. 

Recently, a determination of the gluon density based on single-inclusive charm 



measurements has been performed by the HI Collaboration fllSfl . The results are 
encouraging: the gluon density g(x, fi 2 ) has been determined in the range 7.5 • 10~ 4 < 
x < 4- 10~ 2 , at average scales fi 2 = 25 to 50 GeV 2 , and the result is in good agreement 
with those obtained from the analysis of scaling violations of F2. 

The possibility of performing experiments at HERA with longitudinally polarized 
beams has also been considered | 19| . At leading order in QCD, the heavy-flavour 



production cross section in polarized ep collisions is proportional to the polarized 
gluon density in the proton, Ag. Therefore, data on charm production could in 
principle be used to measure Ag directly, as in the unpolarized case. This possibility 
was first suggested in refs. |2(J and reconsidered in ref. J2l| and in ref. [p2fl ; in this 
last paper the impact of the hadronic component has been estimated. Finally, a full 
NLO computation of heavy quark photoproduction with polarized beams has been 
presented in ref. [p3fl . 

The result is that total cross section asymmetries for the point-like component 
are quite small in absolute value, and can be measured only if very high experimental 
efficiencies and luminosities are achieved; furthermore, the total cross section asym- 
metry is sizeably contaminated by the hadronic process [^2 |. The situation improves 
when considering more exclusive quantities; in ref. |H] it was found that at moder- 
ate p T values the asymmetry for the point-like component can be rather large, well 
above the minimum observable value (in this region, the experimental efficiency is 



sizeable |T5|]). In ref. |22j it was argued that the hadronic component should have a 
negligible impact in this case. 

It can be concluded that heavy-quark photo- and electroproduction studies at 
HERA are still in an early stage. Much more statistics will become available in the 
near future, allowing more accurate tests of QCD. The agreement of present data on 
charm total cross sections with the NLO QCD calculations are satisfactory for choices 
of parameters consistent with low-energy fixed-target results. Transverse momentum 



6 



and rapidity distributions have also been presented, and a general agreement with the 
theory was found. Some discrepancies, as found for example in the rapidity distri- 
butions, are not statistically compelling as yet. With higher statistics and improved 
efficiencies, measurements of the gluon density in the proton will greatly improve. 
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